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Abstract. In its final year of inspiral, a stellar mass (1 — IOMq) body orbits a 
massive (10^ — IO'^Mq) compact object about 10^ times, spiralling from several 
Schwarzschild radii to the last stable orbit. These orbits are deep in the massive 
object's strong field, so the gravitational waves that they produce probe the strong 
field nature of the object's spacetime. Measuring these waves can, in principle, 
be used to "map" this spacetime, allowing observers to test whether the object 
is a black hole or something more exotic. Such measurements will require a good 
theoretical understanding of wave generation during inspiral. In this article, I 
discuss the major theoretical challenges standing in the way of building such 
maps from gravitational-wave observations, as well as recent progress in producing 
extreme mass ratio inspirals and waveforms. 
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Chandrasekhar has described black holes as "the most perfect macroscopic objects 
there are in the universe" . This description refers to their simplicity, depending (in 
astrophysical contexts) solely on the mass and spin of the hole. This dependence, in 
turn, follows from the black hole uniqueness theorems ||; ||, |[ |[ ^, which guarantee 
that all of the other "hairs" will radiate away during the hole's formation. If general 
relativity correctly describes gravity, then the massive compact objects at the centers 
of most galaxies are probably described exactly by the Kerr rotating black hole 
metric. In principle, this description can be tested using LISA Q: gravitational- 
wave observations of "small" (1 — 10 Mq) bodies spiralling into "large" (10^ — 10^ Mq) 
compact objects can be used to map the spacetime of the large compact object, testing 
whether it is a Kerr black hole or some other exotic object. In this article, I discuss 
the theoretical challenges that must be met before such maps can be constructed, as 
well as recent progress. 

Fintan Ryan j|] first showed that a body's spacetime can be mapped with 
gravitational waves. The spacetime of a massive object arises from its multipole 
moment structure. These multipole moments come in two varieties, mass (M;) and 
current (Si). If p{r) is the mass density at position r, and v{r) is the fluid velocity at 
f, then these multipoles are roughly 



M; ~ y dVr'p(r), Si ~ / dVr'-i [fx w(f)p(f)] (1) 
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For a black hole, the mass and current multipole moments take a far simpler form: 

Mf^ + iSf^ = M{iaf (2) 

Because these moments directly determine the spacetime, the orbits of and radiation 
emitted by a smallQ body moving in this spacetime are strongly influenced by the 
massive object's multipoles. Measuring the radiation allows one (in principle) to 
measure the multipolar structure of the massive object. 

Much work remains before it will be possible to measure an exotic body's 
multipole moments in practice. In particular, techniques must be developed to solve 
the wave equation for gravitational radiation from orbits of exotic objects. This is 
not a terribly difficult matter for black hole spacetimes — following Teukolsky's 1972 
discovery ||^ that the wave equation for radiation propagating in Kerr spacetimes is 
separable, an array of calculational technology has been developed for studying the 
generation and propagation of such radiation (see, e.g., |l0[| for review and discussion). 
Very little work has been done for radiation generated in and propagating through the 
spacetime of more exotic objects: Ryan has examined the scalar waves produced by 
highly constrained orbits (circular and equatorial) of objects with arbitrary multipole 
moments. We are a long way from understanding the waveforms created by orbits 
of non-black hole objects. Such an understanding will be needed — at least to some 
degree — in order to probe the multipole character of massive compact bodies. 

For now, we simplify the problem by focusing upon waves generated when the 
massive object is a black hole. Of greatest interest is understanding the waves emitted 
as a small, spinning body spirals into a massive Kerr black hole. Neglecting radiation 
reaction, the motion of such a body is governed by the Papapetrou equations fl^ : 

Dv'^ 1 
Dt " 2^""^'' ^ 

= p^v^^p'^v'^ (3) 

Ut 

In these equations, D/Dt denotes a covariant derivative along the small body's 
trajectory, W^p^ is the Riemann curvature tensor of the background spacetime, S^'^ is 
a tensor related to the spin of the body, = dz'^/dr where z'^{t) is the coordinate 
worldline of the body, and is a generalization of the body's 4-momentum that 
incorporates spin. These equations show that as the small body orbits, its spin 
couples to the curvature of the massive black hole. In the weak field, this coupling 
is known to lead to precessional effects which modulate the phase and amplitude 
of the gravitational waveform In the strong field, the coupling might become 

extremely important. It is likely that, at least for certain parameter values, the orbital 
character will become chaotic. Janna Levin Jl4] has shown that integrating the second 
post-Newtonian equivalents of Equation (^^leads to chaotic motion: orbits with 
similar initial conditions evolve to vastly different configurations. This could make 
data analysis extremely difficult — matched filtering, for example, would require an 
enormous number of templates in order to cover all possible inspirals. When radiation 
reaction is included, the effects of chaos in Levin's work become less extreme: she 
focused upon sources of interest to ground-based detectors, and the number of orbits 
visible to those detectors is not very large. With LISA, though, we expect to see 
around 10^ orbits. Chaotic evolution could lead to a dramatic divergence of outcomes 
from similar initial conditions, rendering data analysis practically impossible. Detailed 
studies of these orbits with parameters relevant to LISA are urgently needed. 

X "Small" means that the orbiting body does not significantly change the spacetime. 



Gravitational waves from extreme mass ratio inspirals 



3 




Figure 1. A slice of the orbital torus. The dark circle is the central Kerr black 
hole; the grey region is the volume which is filled by the orbiting body. The full 
torus is given by rotating this figure about the spin axis (indicated by the arrow). 



Since spin opens an as-yet-poorly-understood can of worms, we will ignore it for 
now, treating the small body as a point perturbation to the Kerr spacetime. At zeroth 
order, this small body moves on a Kerr geodesic. For such orbits, a relatively mature 
radiation reaction formalism (see, e.g., and references therein) has been developed 
that finds the first order radiative corrections to this motion, allowing one to compute 
the trajectory that the body follows as it spirals into the black hole, and the waveforms 
that it generates. This formalism uses the zeroth order geodesic motion as a source 
for the first order corrections. As such, it assumes that the inspiral is adiabatic: the 
timescale for radiation reaction to change the orbit's characteristics is much smaller 
than the orbit's dynamical timescale. 

For the most interesting orbits — eccentric, inclined orbits in the strong field of 
rotating black holes — it is not clear at the present time if this adiabatic approximation 
is reasonable. It is somewhat difficult to define the orbit's dynamical timescale in this 
interesting case. The issue is that there are three physically meaningful timescales: the 
time for the body to cover < </> < 27r; the time for the body to move from rmax 
to r,„in and back; and the time Tq for the body to move from its highest latitude 0min 
to its lowest 6'max and back. These timescales generically are rather different. When 
the orbit is constrained (eccentric but equatorial, or inclined but circular) only two 
of these timescales are non-zero. One can analyze the orbit in a frame that rotates 
at frequency O = cancelling out the <j} motion. It is then simple to define 

the orbit's dynamical timescale: it is (for eccentric, equatorial orbits) or Tg (for 
inclined, circular orbits) (see for more detail). 

When the motion is not constrained, this trick does not work. To understand 
the behavior of the orbit for the general case, consider Figure 0. This figure shows 
a 2-dimensional slice of the volume that is filled by a small body as it orbits the 
black hole: as it moves from r-a^in to rmax it simultaneously moves between ^min and 
^max- Rotating this figure about the black hole's spin axis, we see that the orbit 
ergodically fills a torus in the spacetime near the hole's horizon. For adiabaticity 
to be a good approximation, the orbiting body must come "close" to every point 
in this volume. (The meaning of "close" is of course rather ambiguous. "Close 
enough" depends upon the accuracy that one requires, which in turn depends upon 
how well one needs to know the phase of the waveform.) An important problem is to 
determine how long it takes for an orbiting body to come "close enough" to all points 
in this volume for parameters that are of interest to LISA observations (strong field. 
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large black hole spin, high eccentricity, and arbitary inclination angle). If this time 
turns out to be larger than the radiation reaction timescale then further studies of 
these waveforms will require radiation reaction forces that do not rely on adiabaticity 
@ 0, [9[ |2l|, H. The phase evolution of waveforms that one constructs in 
these cases may have a strong dependence on the orbiting body's initial conditions — 
they may be "effectively chaotic" in the words of Schutz . 

Even if it turns out that the issue of adiabaticity is not a serious problem, 
there is still an important challenge that must be faced in order to evolve generic 
Kerr orbits. Such orbits are characterized by three conserved quantities: the energy 
E, the z-component of angular momentum L^, and the Carter constant Q. For 
Schwarzschild black holes, Q is just L^+L^ , the other angular momentum components. 
Interpretation of the Carter constant becomes less clear as a increases (the geometry 
becomes oblate, confusing the meaning of x, y, and z, and frame dragging entangles the 
t and (j) coordinates), but it is useful to regard it essentially as the "rest" of the small 
body's angular momentum. At present, the most mature computational formalisms 
(such as that described in ||l^) cannot evolve Q. These formalisms work by a method 
of "flux-balancing" : from the flux of gravitational waves going to infinity and down 
the hole's event horizon, one can easily deduce the change in E and because of 
gravitational-wave emission. One cannot easily deduce the change in Q, except in 
special cases (for eccentric equatorial orbits, Q = at all times; for inclined, circular 



orbits theorems which prove that the orbit adiabatically remains circular 24, p5| 



allow one to express dQ/dt as a function of dE/dt and dL^/dt). Although clever 
methods may make it possible to evolve Q just by examining the radiation flux at 
infinity and at the horizon (see Wolfgang Tichy's contribution to these proceedings), 
it may turn out that a local radiation reaction force will be needed. 

At this point, the current unsolved or poorly understood issues have narrowed 
the class of sources that are well understood rather severely. Current computational 
technology is limited to understanding the adiabatic evolution of spinless bodies on 
constrained orbits — either eccentric equatorial orbits or inclined circular orbits. For 
the remainder of this article, I will focus on circular inclined orbits, as described in 
p5| ; Daniel Kennefick and Kostas Glampedakis are developing an analysis of eccentric 
equatorial orbits using a similar formalism. 

Using the formalism and code described in jl^, I have studied the inspirals and 
associated gravitational waveforms for a large number of strong-field initial conditions. 
The results for a — 0.998Af are shown in the left-hand panel of Figure ||. This figure 
shows the inclination angle t of the orbit as the small body spirals from r = 4M to 
the last stable orbit (LSO); it also shows the number of days that pass along this 
sequence. In all cases the trajectory is nearly flat — the inclination does not change 
very much as the small body inspirals. The inclination evolution is even flatter at 
smaller values of a Curt Cutler has suggested that this might be used as the 

basis of an approximate scheme for evolving the Carter constant for generic orbits: if 
it is true in general that the inclination angle does not change very strongly, then it 
might be a reasonable approximation to set the change to zero. This condition would 
constrain the evolution of Q. Cutler's approximation may be useful for developing 
approximate waveforms for the study of data analysis tools. 

Turn now to the right-hand panel of Figure ^. This panel is identical to the 
left-hand panel except that the flux of radiation down the massive black hole's event 
horizon has been ignored in constructing the inspiral trajectory. Although the shape 
does not change very much, the time it takes to inspiral is significantly smaller, the 



Gravitational waves from extreme mass ratio inspirals 



5 



t = 35.8 7963 days 
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Figure 2. Inspiral from r = AM to the last stable orbit (indicated by the dotted 
line). The massive black hole has M = 10^ Mq and a = 0.998M; the small 
body has 11 = 1 Mq . The left-hand panel shows the inspiral including both flux 
to infinity and flux down the horizon; the right-hand panel includes only flux to 
infinity. Notice that the horizon flux typically slows the inspiral. 



horizon flux slows the inspiral by several weeks at low inclination angle. This is a 
significant effect — change of the inspiral time by such a large amount should be 
easily measurable. At first glance, it is also rather counterintuitive: one expects the 
hole's event horizon to be a sink of energy, so the inspiral would speed up as radiation 
flows into the hole. This simple picture is wrong when the hole rotates. A more 
accurate picture can be developed by considering the tidal coupling of the black hole 
to the inspiralling body. The tidal field of the small body will distort the hole, raising 
"bulges" in the event horizon These bulges exert a torque back on the small 

body. When the hole is rapidly spinning, the bulges are dragged ahead of the orbiting 
body so that this torque tends to increase the orbiting particle's energy. This partially 
offsets the energy that is lost from radiation to infinity, slowing the inspiral. 

One of the major goals of this analysis is to produce gravitational waveforms. 
Using the inspiral trajectories shown in the left-hand panel of Figure ^, I have 
developed the associated gravitational waveforms: 

h^it) - i/.x (t) = ^ E E E SU^) eW™—*)] (4) 

i=2 m=-i fc=-oo ™« 

Here, /i is the mass of the small body, D is the luminosity distance to the source, 
[r{t),t{t)] is the inspiral trajectory, Zimk[r,L] is a complex amplitude computed with 
the Teukolsky equation, ujmk = milfi, -\- fcfie, Sim is a spin- weighted spheroidal 
harmonic, and if) is the angular position of the observer relative to the spin axis. 
Perhaps the most effective demonstration of the characteristics of these waveforms is 
given by converting the functions h+{t) and hx{t) into sounds; the reader is invited 
to visit the URL given in and listen to the sounds available there. Some of the 
features one can hear in these waveforms are rather surprising. For instance, in several 
cases, the wave chirps down as well as chirps up: a portion of the sound has decreasing 
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frequency. Despite the many simplifications that were imposed in order to compute 
these waves, they have a rather complex and ornate character. 

As the challenges discussed earlier are surmounted we will be able to develop 
waveforms that incorporate even more structure and complexity. The surprising 
features that were found for simple circular orbits will doubtless be joined by more 
surprises, increasing the complexity of the waveforms. This complexity will make it 
difficult to detect and analyse the waves in LISA data, but is indicative of how much 
can be learned from their observation. 
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